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INTRODUCTION 


It  is  well  known  that  plasticity  plays  an  important  role  in  crack 
growth  in  metals.  Due  to  the  pronounced  plastic  behavior  in  many  of 
the  matrix  metals  (e.g.  aluminum  or  titanium)  used  in  metal-matrix 
composites  the  same  is  expected  to  be  true  for  these  new  materials  also. 
Consequently,  establishing  the  relation  between  the  matrix  plasticity 
and  the  overall  fatigue  properties  of  the  composite  will  aid  conside’'ably 
in  understanding  the  crack  growth  process. 

One  of  the  most  important  plasticity  effects  in  fatigue  is  the 
retardation  of  the  crack  growth  due  to  overloads.  There  have  been 
extensive  studies  on  the  overload  effects  in  homogeneous  metals,  but  no 
research  has  been  done  on  this  aspect  of  the  fatigue  behavior  in  metal- 
matrix  composites.  The  present  research  is  part  of  a  systematic  study 
to  obtain  such  information. 

No  meaningful  study  of  retardation  effects  can  be  pursued  without 
first  performing  an  in-depth  analysis  of  the  elastic-plastic  behavior 
of  the  material.  Consequently,  considerable  effort  was  spent  in  this 
stage  and  the  present  report  summarizes  the  results. 

Chapter  1  describes  some  of  the  preliminary  tests  performed  on  the 
B/Al  composite.  It  was  concluded  that  this  material  appeared  to  oehave 
like  an  orthotropic  elastic-plastic  material.  Chapter  2  then  proceeds 
to  a  proper  formulat  m  of  anisotropic  elastoplasticlty  plus  a  scheme 
for  determining  the  relevant  material  properties.  Chapter  3  shows  the 
results  of  the  experimental  characterization.  Finally,  Chapter  4  considers 
the  problem  of  deterrr'nlng  experimentally  the  stress  and  strain  fields 
surrounding  the  crack  tip  during  fatigue  loading. 
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The  studies  reported  here  form  the  foundation  for  the  later  experi¬ 
mental  and  numerical  studies  of  the  effects  of  plasticity  on  the  fatigue 
crack  propagation  in  B/Al  metal-matrix  composites. 
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Chapter  1  PRELIMINARY  TESTS 


In  order  to  make  a  decision  as  to  how  to  proceed  with  the  character¬ 
ization  of  a  metal-matrix  material  a  working  knowledge  of  its  basic 
behavior  had  to  be  established.  First,  tensile  tests  were  done  to 
determine  principal  direction  m.  terial  properties.  Additionally,  some 
fatigue,  fracture,  and  overload  tests  were  also  conducted  to  explore  its 
major  crack  growth  behavior. 

MATERIAL  PROPERTIES 

The  two  plates  of  boron /aluminum  received  from  AMERCOM  were  nominally 
12"  X  24".  These  plates  consisted  of  11  plies  of  5.6  mil  boron  fibers 
in  a  matrix  of  6061  aluminum.  The  ply  layup  was  unidirectional  for  both 
plates,  the  fibers  running  lengthwise.  The  nominal  thickness  of  the 
plates  was  .077",  and  the  fiber  volume  fraction  was  47.5%.  Both 
plates  showed  signs  of  being  slightly  warped  across  the  direction 
perpendicular  to  the  fibers. 

The  specimens  were  cut  from  the  panel  using  a  6”  diameter  diamond 
wheel  .04"  in  thickness,  turning  at  1335  RPM.  Water  cooling  was  employed 
to  protect  the  material  from  frictional  heating.  The  cutting  machine 
had  an  automatic  feed  table  to  which  the  panel  was  clamped.  The  feed 
rate  was  approximately  5  inches  per  second  back  and  forth,  and  the  blade 
was  slowly  lowered  through  the  thickness  of  the  panel  at  a  rate  of 
approximately  .001"  per  eight  passes  of  the  table.  For  cutting  notches 
in  the  fatigue  and  fracture  specimens,  the  same  machine  was  employed 
with  the  edge  of  the  specimen  held  perpendicular  to  the  blade  and  with 
the  table  stationary. 
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Figure  1  shows  the  specimen  cutting  plan  for  the  first  of  the  two 
panels.  Table  1  is  a  list  of  all  specimens  cut  from  panel  1,  giving 
their  dimensions  and  uses.  Five  specimens  of  nominal  dimensions  9"  x  1" 

X  .077"  were  tested  using  the  setup  shown  in  Figure  2  to  give  the  ortho- 
tropic  material  properties.  Aluminum  tabs  were  used  on  all  specimens 
to  protect  the  fibers  from  being  crushed  by  the  waffle  pattern  of  the 
grips.  The  tabs  were  applied  with  cyanoacrylate  glue  (commonly  referred 
to  as  "super-glue") .  The  strength  of  the  glue  is  superior  to  that  of 
epoxy,  and  much  more  convenient  to  use  because  it  sets  up  in  seconds 
rather  than  hours. 

Clip  Gage 

The  strains  were  taken  from  the  specimens  using  a  clip  gage  manu¬ 
factured  by  Measurements  Technology,  Inc.,  and  it  allowed  both  longitu¬ 
dinal  and  transverse  strains  to  be  recorded  simultaneously.  Four  arms 
ending  in  pins  attach  the  gage  to  the  specimen  as  shown  in  Figure  3. 

Two  of  the  arms  bend  freely  as  the  specimen  elongates,  and  strain  gages 
on  these  arms  yield  the  longitudinal  strain.  As  the  specimen  width  is 
changed  by  the  Poisson's  contraction,  the  transverse  strain  is  picked 
up  by  the  gages  on  the  spring  which  holds  the  pins  onto  the  specimen. 

During  testing,  both  longitudinal  arms  of  the  clip  gage  were  wired  into 
a  Wheatstone  bridge  to  complete  a  full-bridge.  The  single  transverse 
arm  of  the  gage  was  wired  into  a  half -bridge.  In  all  cases  where  the 
clip  gage  was  used,  a  10  volt  bridge  excitation  was  maintained.  Averaging 
both  longitudinal  arm  sensitivities,  the  bridge  output  voltage  is  con¬ 
verted  to  specimen  strain  by  the  factor  1  rav  =  1155  ye.  For  the  transverse 
direction  the  strain  conversion  factor  is  1  mv  =  3650  ye.  (Note  the  lack  of 
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sensitivity  in  the  transverse  direction.)  Figure  3  shows  a  comparison 
of  both  strain  gage  and  clip  gage  results.  The  data  for  this  particular 
comparison  comes  from  the  test  conducted  on  B/Al  specimen  5  (10®  off- 
axis  specimen) .  The  resistance  gages  and  the  clip  gage  both  show  essen¬ 
tially  the  same  results.  The  clip  gage  has  the  advantage  of  measuring 
the  strain  of  the  middle  surface  and  consequently  will  not  pick  up  out- 
of-plane  bending.  This  is  significant  since  the  panels  exhibited  some 
initial  warping.  Further,  the  clip  gage  is  more  reliable  for  obtaining 
the  large  strains  after  yielding.  The  particular  clip  gage  used  was 
calibrated  at  the  factory  and  was  shown  to  give  linear  output  up  to  a 
strain  of  ,025  in/in.  The  compliance  of  the  clip  gage  was  not  a  problem 
due  to  the  relative  stiffness  of  the  material. 

Combining  the  load  voltage  output  from  the  MTS  810  control  panel 
and  the  Wheatstone  bridge  outputs,  X-Y  plotters  were  used  to  record 
the  load-strain  histories.  The  following  sections  provide  details  of 
the  five  uniaxial  specimen  tests  and  the  results  obtained.  A  summary 
of  the  orthotropic  material  constants  obtained  from  these  tests  is 
shown  in  Table  2.  A  plot  of  the  typical  principal  material  direction 
stress-strain  curves  obtained  from  these  tests  is  shown  in  Figure  4. 

Note  the  ductility  in  the  transverse  direction  and  nearly  brittle  be¬ 
havior  in  the  longitudinal  direction. 

0°  Tests 

Specimen  6  was  loaded  monotonlcally  to  failure,  and  only  longitudinal 
strain  data  was  taken.  The  material  In  this  test  displayed  nearly 
linear  behavior  up  to  fracture.  (Although  the  tests  on  specimens  6  and 
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7  were  conducted  with  the  MTS  on  load-control,  later  tests  run  under 


stroke-control  gave  identical  results.)  The  ultimate  stress  was  222  ksi. 
The  initial  elastic  modulus  was  found  to  be  30.3  msi. 

Specimen  7  was  tested  on  two  occasions;  the  first  test  was  to  obtain 
the  elastic  material  constants  (highest  stress  level  of  this  test  was 
40  ksi)  and  the  second  test  conducted  was  a  series  of  loading  and 
unloar’ings  to  observe  the  cyclic  hysteresis  behavior.  Up  to  a  stress 
of  20  ksi  the  average  slope  of  the  stress-strain  curve  yielded  a  value 
of  32.6  msi  for  Young's  modulus.  At  a  stress  level  of  40  ksi  this  value 
dropped  to  29.5  msi.  In  the  second  test  the  specimen  was  loaded  in 
cycles  of  increasing  maximum  stress  in  order  to  observe  whether  there 
was  any  residual  strain  when  unloaded.  Despite  the  fact  that  the  unloading 
curves  never  followed  the  loading  curves,  the  strain  always  returned  to 
zero;  no  residual  strain  was  observed.  In  the  first  loading  cycle, 
longitudinal  and  transverse  strains  yielded  a  value  of  .29  for  Poisson’s 
ratio  ('^2^2^*  each  loop,  an  average  elastic  modulus  (E^  was  cal¬ 

culated  using  the  two  reversal  points  for  the  slope.  The  test  and 
results  are  summarized  below: 

9  cycles  at  maximum  stress  of  160  ksi 
Average  modulus  E  =  29.4  msi 

5  cycles  at  maximum  stress  of  190  ksi 
Average  modulus  E  =  28.9 

1  cycle  at  maximum  stress  of  196  ksi 
Average  modulus  E  =  28.5  msi 

Cycle  up  to  failure  at  stress  of  204  ksi 
Figure  5  shows  a  reproduction  of  the  last  eight  cycles  before  failure 
of  specimen  7  taken  from  the  X-Y  plotter.  Each  cycle  was  spaced  apart 
for  clarity. 
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90“  Tests 


Specimen  9  was  loaded  monotonically  to  failure.  Extensive  yielding 
was  observed,  with  a  failure  strain  of  .0132  in/in  and  an  ultimate  load 
of  27.0  ksi.  Yielding  was  first  observed  at  a  stress  of  14.5  ksi.  Using 
a  .05%  yield  criterion,  the  yield  stress  was  16.8  ksi.  The  failure 
stress  of  this  specimen  was  26.5  ksi.  A  value  of  V2j  =  -37  was  calculated 
from  the  transverse  strain  data  in  the  linear  region. 

Specimen  8  was  loaded  to  failure  but  with  intermittent  unloadings 
to  observe  the  residual  strains  associated  with  the  yielding.  Figure  6 
is  a  reproduction  of  the  plotted  data.  From  the  elastic  loading  slope 
the  value  of  Young's  modulus  was  calculated  to  be  20.3  msi  while  a  value 
of  .21  was  calculated  for  ^2^.  The  onset  of  yielding  occurred  earlier 
than  it  did  in  specimen  9  at  a  stress  of  10.3  ksi.  According  to  the 
.05%  yield  criterion,  the  yield  stress  was  15.6  ksi.  Similar  to  the 
cyclic  test  of  specimen  7,  the  unloading  paths  were  different  than  the 
loading  curves.  This  time,  of  course,  residual  strains  were  observed. 

For  each  unloading  cycle  a  value  of  average  modulus  was  calculated  (again 
using  the  two  reversal  points)  and  these  values  were  plotted  against 
the  plastic  strain  at  which  they  occurred.  The  plot  of  this  data  can 
be  seen  in  Figure  7.  Note  that  the  material  appears  to  soften  with 
progressing  strain.  This  softening  is  probably  due  to  fiber/matrix 
debonding. 

10°  Test 

The  purpose  of  doing  an  off-axis  test  was  to  obtain  a  value  of  the 
elastic  shear  modulus  and  an  approximate  value  for  the  ultimate  shear 
stress.  The  idea  for  a  10“  specimen  was  obtained  from  Reference  [1], 
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This  paper  recommends  that  the  10“  off-axis  tensile  specimen  be  considered 
as  a  possible  standard  for  characterization  of  inplane  shear  of 
unidirectional  fiber  composites.  The  method  requires  strain  data  from 
a  gage  rosette  to  calculate  the  inplane  shear  strain,  while  the 
Inplane  shear  stress  is  found  from  the  stress  transformation 
equation.  The  slope  of  the  resulting  shear  stress-strain 
curve  yields  the  shear  modulus  necessary  equations  for  this 

reduction  are 


e  =  e, 

X  1 

e^  =  0“  gage 

e^  =  (2e2  +  2e^  -  e^)/3 

=  120“  gage 

V  '  2(^3  - 

e^  =  240“  gage 

e,-  =  (e  -  e  )sin20  +  e  cos26 

12  y  X  xy 

and  0  =  10“ 

0.^  -ho  sin26 

12  X 

(1.1) 

The  value  of  the  elastic  shear  modulus  obtained  by  this  method  from 

specimen  5  was  7.55  msi.  The  shear  modulus  was  also  calculated  from 

the  aquation  of  the  apparent  elastic  modulus  for  an  orthotropic  lamina 

stressed  in  a  nonprincipal  material  direction.  The  equation  for  the 

apparent  elastic  modulus,  which  is  good  only  for  uniaxial  stress,  is 

1/Eg  =  cos^0/Ej^  +  sin^9/E2  + 

(1/Gj^2  '  2Vj2/Ej)sin^ecos^0  (1.2) 

Solving  for  the  value  of  the  elastic  shear  modulus  obtained  by 

this  method  was  7.65  msi,  which  compares  favorably  with  the  value  from 

the  slope  of  the  stress-strain  curve.  A  parametric  study  of  the  above 

equation  reveals,  however,  that  the  value  of  Gj2  calculated  is  very 

sensitive  to  the  values  of  E„  and  E,  used. 

y  1 

The  value  of  the  intralaminar  shear  stress  0^2  was  calculated  to 
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be  19. A  ksi  at  specimen  failure.  Figure  4  does  not  show  the  failure  point 
of  the  shear  curve  because  it  is  off  the  graph  at  approximately  .083  in/in 
shear  strain. 

FATIGUE  TESTS 

All  fatigue  tests  were  performed  using  the  setup  shown  in  Figure  8. 
The  microscope  used  was  made  by  the  Precision  Tool  and  Instrument  Co. 

Ltd.  and  gave  an  approximate  x7  magnification.  Cross  hairs  in  the  eye¬ 
piece  were  used  to  locate  the  crack  tip  and  the  relative  movement  of 
the  crack  was  read  off  the  vernier  scales  of  the  microscope.  Each  test 
was  run  with  a  constant  amplitude  sinusoidal  load  cycle,  tension-tension 
(Omin  -  0),  at  a  frequency  of  10  Hz  at  room  temperature. 

To  aid  visualization  of  the  crack  tip,  Spotcheck  dye  penetrant 
made  by  Magnaflux  was  applied,  along  with  a  developer.  Tests  on  regular 
aluminum  specimens  were  conducted  with  and  without  the  Spotcheck  penetrant 
to  determine  whether  it  had  any  effect  on  the  crack  growth  rate.  In 
the  low  cycle  fatigue  range  typical  of  these  tests  there  was  no  indication 
that  the  penetrant  had  any  effect  on  the  crack  growth  rate. 

Crack  measurements  were  made  with  the  machine  stopped.  These 
measurements  were  taken  only  on  one  side  of  the  specimen  which  had  been 
polished  prior  to  cycling  with  600  grit  paper  and  finished  with  tripoli 
buffing  compound. 

Note  that  the  notches  cut  into  the  specimens  to  initiate  the  fatigue 
cracks  were  at  least  .04"  wide  (the  width  of  the  diamond  blade).  The 
notch  roots  were  not  quite  square  but  were  somewhat  rounded. 
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90“  Test 


For  specimen  10,  the  crack  initiated  from  the  notch  root  and  pro¬ 
pagated  parallel  to  the  fibers  (perpendicular  to  the  loading  direction). 

The  results  are  similar  in  nature  to  what  would  be  expected  from  a  crack 
growing  transverse  to  the  load  in  a  homogeneous  material. 

0°  Tests 

For  the  specimens  with  a  notch  cut  perpendicular  to  the  fibers  it 
was  found  to  be  impossible  to  initiate  a  transverse  fatigue  crack. 

Specimen  1  was  given  a  .166"  notch  and  cycled  at  a  maximum  stress  of 
22  ksl,  which  was  not  sufficient  to  grow  a  crack  through  the  fibers. 

Figure  9  shows  in  schematic  form  how  the  cracks  propagated  along  the 
fibers  in  a  mode  II  fashion.  The  notch  was  then  cut  deeper  into  the 
specimen  (from  .166"  to  .315")  and  again  a  transverse  crack  could  not 
be  initiated.  Cutting  the  notch  back  even  further  (from  .315"  to  .74") 
and  cycling  at  a  maximum  stress  of  15.2  ksi,  the  cracks  were  again  ob¬ 
served  to  only  grow  parallel  to  the  fibers.  From  this  point  on,  the 
nominal  width  of  3"  for  the  specimens  was  thought  to  be  wasteful  for 
the  0“  specimens.  Specimen  1  itself  was  cut  into  two  pieces  to  perform 
two  tests:  Specimen  la  with  a  .74"  notch  and  the  longitudinal  fatigue 
cracks  still  in  it  was  loaded  to  failure  to  observe  the  failure  mode, 
and  specimen  lb  was  given  a  notch  of  .122".  The  results  of  the  test 
on  specimen  la  are  given  in  the  following  section  of  fracture  tests. 

Specimen  lb  was  cycled  at  a  maximum  stress  of  37.1  ksi  and  the  cracks 
that  initiated  at  the  notch  root  again  grew  parallel  to  the  fibers.  The 
load  was  Increased  at  this  point,  but  even  up  to  the  cyclic  load  of  79.3  ksi 
the  cracks  did  not  turn  but  grew  rapidly  along  the  fibers.  The  notch  was 
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cut  deeper  (from  .122"  to  .281")  and  given  an  initial  cyclic  load  of 
61.7  ksi.  Again  the  cracks  initiated  from  the  notch  tip  and  traveled 
parallel  to  the  fibers. 

In  addition  to  providing  crack  growth  data  specimen  15  was  used  to 
provide  Information  on  the  relative  load  levels  required  to  propagate 
cracks  in  mode  II.  For  the  specimen  geometry  given  in  Table  1,  the 
fatigue  threshold  was  found  to  be  approximately  8.5  ksi,  where  da/dN  was 
calculated  to  be  nearly  300  x  10  in/cycle. 

Off-Axis  Fatigue  Test 

To  confirm  that  the  fatigue  cracks  In  boron/aluminum  will  always 
propagate  parallel  to  the  fibers,  a  15®  off-axis  specimen  (Specimen  2)  was 
machined  so  that  the  notch  was  perpendicular  to  the  load  (and  75®  to 
the  fibers).  Once  again  the  crack  grew  parallel  to  the  fibers  when  a 
cyclic  stress  of  10.8  ksi  was  applied. 

FRACTURE  TESTS 

The  fracture  loads  for  the  material  at  different  fiber  orientations 
were  required  in  order  to  know  how  much  of  an  overload  could  be  applied 
to  the  fatigue  specimens. 

The  fracture  test  conducted  on  specimen  la  showed  that  a  0*  specimen 
with  pre-grown  cracks  (parallel  to  the  fibers)  still  has  considerable 
strength.  A  load  of  8230  lbs.  was  required  to  break  specimen  la.  This 
corresponds  to  a  remote  stress  of  55.7  ksi,  or  a  net  section  stress  of 
90.6  ksi.  The  fatigue  cracks  In  the  specimen  were  observed  to  extend 
only  a  small  amount  before  the  specimen  ruptured.  The  fracture  surface 
ran  perpendicular  to  the  fibers.  Since  no  cracks  could  be  Initiated 
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through  the  fibers  in  the  0®  direction,  all  that  could  be  done  was  a 
fracture  test  of  a  notched  specimen  without  fatigue  cracks.  This  was 
done  on  specimen  14.  The  notch  radius  was  approximately  .025  inch. 

The  fracture  load  was  8775  lbs.,  corresponding  to  a  remote  stress  of  74  ksi. 

Fracture  occurred  in  specimen  13  while  attempting  to  apply  an  over¬ 
load  of  2.0  times  the  maximum  cyclic  stress.  The  crack  length  had 
grown  to  .301  inch  when  the  overload  was  attempted.  The  maximum  load 
achieved  was  4150  lbs.  (17.56  ksi  remote  stress). 

DISCUSSION 

Concerning  the  cutting  of  the  materials,  because  of  the  relative 
brittleness  of  the  boron  fibers  it  is  possible  that  the  diamond  saw 
may  have  caused  local  damage  to  the  fibers.  Since  the  cutting  was  water 
cooled,  heating  of  the  specimen  was  not  a  problem.  An  altogether  dif¬ 
ferent  method  of  cutting  the  material  with  electric  discharge  machining 
is  introduced  in  Chapter  4. 

Since  the  panels  delivered  had  rough  edges,  it  was  not  possible  to 
determine  the  exact  fiber  orientation  in  the  plate  when  laying  out 
the  cutting  lines.  It  was  fortunate  that  in  the  manufacturing  process 
the  fibers  caused  a  ridged  pattern  to  be  formed  on  the  surface  of  the 
plates  by  which  the  fiber  orientation  could  be  found.  To  check  the 
accuracy  of  the  cutting,  the  following  analysis  was  performed;  Close-up 
photographs  of  the  edges  of  the  specimens  exposed  by  the  cut  were  taken. 

The  pictures  showed  a  light  and  dark  contrast  between  the  matrix  and  the 
fibers.  The  fibers,  which  appeared  as  dark  patches,  formed  ellipses 
when  cut  at  an  angle.  By  taking  measurements  from  the  photographs  of 
the  length-to-width  ratios  of  the  exposed  fiber  face,  the  angle  of  the 


cut  could  be  determined  as  shown  in  Figure  10.  Note  that  the  scale  of 
the  photograph  was  not  important  since  the  dimensions  were  ratioed. 

Also  note  that  this  analysis  is  only  accurate  for  shallow  angles  (i.e. 
for  specimens  with  fiber  orientations  about  10“) .  The  results  obtained 
indicated  that  the  cutting  was  accurate  to  within  a  one  degree  tolerance. 

During  the  course  of  testing,  certain  advantages  and  disadvantages 
were  observed  when  using  the  clip  gage  in  comparison  with  strain  gages. 
These  are  summarized  below 

Advantages  Disadvantages 

•  measure  strain  of  •  lack  of  sufficient 

neutral  surface  sensitivity  in  transverse 

direction 


•  two  longitudinal  arms 
average  Inplane 
bending 

•  ease  and  speed  of 
gaging  specimens 

•  reliable  for  measuring 
high  strain  values 


•  difficult  to  mount  if 
specimen  is  too  thin 

•  can  only  measure  longi¬ 
tudinal  and  transverse 
strains 


Concerning  the  off-axis  specimens  used  to  calculate  the  inplane 
shear  characteristics;  The  10"  specimen  was  convenient  in  that  it  did 
not  use  up  an  excess  of  material,  but  it  was  found  that  for  the  metal- 
matrix  composite  analyzed  here  the  angle  of  10"  was  not  optimum  for 
obtaining  the  shear  values.  According  to  Reference  [1],  for  a  10" 
specimen  to  be  an  adequate  means  of  characterizing  the  inplane  shear 
in  a  uniaxial  composite,  the  shear  stress  a^2  n>nst  be  the 
only  coiiponent  of  stress  near  its  critical  value,  and  at  failure  the 
fracture  must  occur  at  the  10"  plane  when  a ^2  reaches  this  critical  value. 
Figure  11  is  a  plot  of  the  normalized  strain  versus  fiber  angle  for  a 
typical  graphite  epoxy  material.  As  can  be  seen,  the  normalized  intralaminar 
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shear  strain  is  a  maximum  at  about  the  10“  orientation,  which  is  a 
desirable  feature  for  this  test  method.  The  orientation  at  which  the 
normalized  shear  strain  reaches  a  maximum  is  a  function  of 
the  orthotropy  of  the  material.  The  paper  cites  results  for  three 
materials  with  ratios  of  31.2,  16.5,  and  3.5;  for  the  latter 

material  the  normalized  shear  strain  reached  its  peak  around  15“ .  For 
the  boron /aluminum  material,  on  the  other  hand,  the  ratio  of  principal 
elastic  moduli  is  1.5,  and  the  peak  is  reached  near  40“  as  can  be  seen 
in  Figure  12.  It  can  be  concluded  that  for  unidirectional  composites 
with  low  orthotropy  as  exhibited  in  these  two  metal-matrix  composites,  the 
10°  off-axis  specimen  is  not  ideal  for  characterization  of  inplane 
shear  characteristics. 

Table  2  includes  values  of  the  material  constants  which  have  been 
reported  in  other  literature  [Ref.  2].  Most  values  compare  favorably, 
although  the  values  of  Poisson's  ratios  do  not  agree  as  well.  Note  that 
for  the  values  found  here,  and  even  for  values  reported  elsewhere,  the 
reciprocal  relation  “  '^21^^2  orthotropic  materials  is  not 

satisfied.  This  Indicates  that  either  the  material  is  not  strictly 
orthotropic  or  the  reported  values  are  not  precise. 

Finally,  the  unloading  tests  performed  confirm  that  this  metal- 
matrix  composite  exhibits  plasticity  in  the  classical  sense.  That  is, 
the  unloading  is  elastic  with  a  slope  nearly  that  of  the  initial  Young's 
modulus.  Further,  the  reloading  curve  is  elastic  up  to  the  previous 
maximum  stress  and  then  after  yielding  work -hardening  occurs  at  about 
the  previous  rate.  This  plasticity,  however,  is  very  different  at  different 
fiber  orientations,  being  nearly  zero  for  the  longitudinal  specimen. 
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Thus,  it  appears  that  the  nonlinear  behavior  of  boron /aluminum  metal- 
matrix  composites  can  be  modeled  in  terms  of  anisotropic  plasticity  with 
work -hardening . 
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Chapter  2  FORMULATION  OF  ORTHOTROPIC  ELASTOPLASTICITY 


The  elastoplastlc  description  of  the  orthotropic  material  behavior  is 
broken  into  two  parts,  the  elastic  modeling  and  the  plastic  modeling. 


ELASTIC  ORTHOTROPIC  FORMULATION 

The  analysis  begins  with  the  generalized  Hooke's  law  for  anisotropic 
materials  with  its  fourth  order  compliance  tensor 

®ij  "  ^ijkl\l  (2.1) 

where  e .  ,  is  the  strain  tensor,  a,,  is  the  stress  tensor  and  S.  .is 
ij  kl  ijkl 

the  compliance.  This  relation  in  contracted  form  is  reduced  to  the 
three  by  three  matrix  equation  given  below  for  the  case  of  an  orthotropic 
material  under  plane  stress 


(2.2) 


The  subscripts  1,2  denote  principal  material  directions.  (See  Reference  [3] 
for  more  details).  We  now  define  the  following  two  new  terms;  the 
effective  strain  e  and  the  effective  stress  a.  These  factors  are  intro¬ 
duced,  along  with  new  s-terras,  by  dividing  the  equations  through  by  $22 
to  yield  the  new  equations  as  shown 


(de/do) 
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where 


de/da  =  S22  =  I/E2 


®11  = 


=  ^2/^12 


®12  ~  “''12^2'^^!  “  “'’21 


(2.3) 


By  using  the  transformation  equations  for  stress  and  strain  for  arbitrarily 
oriented  uniaxial  specimens  (See  Figure  13)  the  following  equation  is 
derived  for  the  effective  elastic  modulus  E 

—  — .  —  4  4 

E  =  do/de  =  [sjjCos  0  +  sin  0 

+  2(Sgg  +  Sj2)sin^8co8^0] (dOg/deg)  (2.4) 

where  do./de.  =  E.  is  the  Yoving's  modulus  for  the  specimen.  Similarly, 

y  o  y 

the  expression  for  Poisson's  ratio  can  be  found  in  terms  of  the  s-values 
as 

2  2  4  4 

^  ^®T0  ^^11  ^  ®  Sj2(8ln  0  +  cos  6) 

0  ”  4  4  2  2 

0  SjjCos  0  +  sin  0  +  2(Sj2  +  ®66^®^  ®  (2.5) 

The  reason  for  defining  the  s-values  will  become  apparent  in  the  following 
section  for  the  plastic  orthotropic  formulation. 

PLASTIC  ORTHOTROPIC  FORMULATION 

The  derivation  of  the  orthotropic  plasticity  relations  begins  with 
the  statement  of  the  flow  rule  for  plasticity  [^]  as 


dePj,j  -  3g/9<T^j  dX* 


(2.6) 


Note  that  the  superscript  'p'  indicates  plasticity.  The  plastic  flow 


function  g  is  assumed  quadratic  in  as 
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de 


P 


ij 


{c‘ 


,  o,  ,  +  C 


ijkl  kl  innij  ran 


.a  }dA* 


=  2C^ 


(2.7) 


This  can  be  written  In  matrix  form  as 


de 


11 


de 


22 


de 


33 


de 


31 


de 


23 


de 


12 


P 


*^11  ^^12 . ‘^16 

*^21  ^^22 . 

‘^31 . 


. *^66 


dX* 


(2.8) 


where  is  symmetric  (i.e.  it  contains  21 
all  of  which  may  be  functions  of  strain) .  If 
orthotropic  (that  is  to  say,  if  the  body  has 
planes  of  material  symmetry)  then  the  matrix 
additional  symmetries.  The  matrix  becomes 


Independent  coefficients 

the  material  is  plastically 

three  mutually  perpendicular 

c^ . .  will  exhibit  certain 
ij 


tcP] 


showing  that  to  describe 
coefficients  are  needed. 


‘^^ll 

^^^12 

cP 
^  13 

0 

0 

0 

cP 

23 

0 

0 

0 

'^3 

^-23 

cP 

33 

0 

0 

0 

0 

0 

0 

,P 

'  44 

0 

0 

0 

0 

0 

0 

0 

0 

_ 

0 

0 

0 

0 

orthotropic  plastic  flow  only  nine 
For  plane  stress  situations  where 


(2.9) 


independent 


033  -  0. 
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*^23  ~  ^  incremental  relations  further  simplify  to 


de' 


11 


de* 


22 


'=‘’l2'’ll  +‘=''22‘’22'> 


de* 


12 


'66  “12 


(2.10) 


Note  that  the  above  equations  are  over-specified  with  the  four  c-values 
and  dX*.  The  dX*  is  retained  so  that  the  formulation  resembles  classical 
plasticity. 

It  is  noted  that  the  change  of  volume  in  the  body  due  to  plasticity 
is  given  by  the  expression  (for  plane  stress) 


dv  =  e^  +  e^  +  e^ 
av  e  11  ^  ®  22  ®  33 


(c^ll  +  c^^2  ^^^12  ^11  ^ 


(2.11) 


The  change  of  volume  will  not  be  considered  further  because  that  would 
involve  the  coefficient  c^i^*  ^^23*  ^^33  occur  in  the  incre¬ 

mental  flow  rule  (2.10). 


Three-Parameter  Formulation 

Let  the  following  definitions  be  introduced; 


r 


11 


^66  =  11 
^^12  "  ^ \1^^  11 


dX  =  cP22dX* 


Then  the  new  Incremental  expressions  are 


(2.12) 
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1  \ 

^ii'^ii  ^12^22) 

>  -  ^ 

^12'^ll  °22  ; 

r®^2' 

1  i 

^6  ^12  J 

dX 

(2.13) 

The  factor  of  proportionality  dX  is  solved  for  from  work  considerations. 
Parenthetically  this  flow  rule  could  have  been  derived  directly  by  taking 
the  flow  function  g  as 

g  =  o  22  12  ^^12°11°22^ 


Effective  Stress  and  Strain 

In  general,  the  plastic  work  increment  is 

dW^  =  a . .de^, . 
ij  Ij 

For  plane  stress  conditions  this  expands  to 

dwP  «  a^jdeP^j  +  a22‘^«^22  ^^12^^^12 


(2.14) 


(2.15) 


Substituting  for  the  incremental  strains  in  terms  of  the  stresses  from 
equation  (2.13),  the  work  Increment  can  be  rewritten  in  the  form 

dVfP  =  +  022  ^^66^^12  ^’^12'^ll'^22^‘^^ 

=  2gdX  (2.16) 


Let  the  effective  stress  and  effective  plastic  strain  be  defined  such 
that 


ad"?  =  dW^  =  2gdX 


(2.17) 


Further,  define  the  effective  stress  as 


30=  r^^Ojj  +  022  ^''66°12  '''  ^^12'^11°22 


2g 


(2.18) 


20 


Therefore 


ad^  =  o^dX 

or  dX  =  I  idP/a)  =  |  (d? /do)  (do/ a)  (2.19) 

The  effective  plastic  strain  consistent  with  the  above  definition  of 
effective  stress  is 

di^  =  -|  odX 

or  [di^]^  =  I  odi^dX  =  I  2gdX^  =  I  o^dX^  (2.20) 

Substituting  into  "o  for  o^^,  022*  and  0^2  in  terms  of  the  plastic  strains 
gives 

[diPl2=|(rii-ri2V'  [dePii^  -f  riideP22^ 

”  ^^12  ^^66^  (2.21) 

Note  that  for  the  special  cases  of  uniaxial  stress  applied  in  the  principal 
material  directions  there  is  no  shear  term  0^2  and  the  expression  for 
the  effective  strain  and  stress  reduces  to 

Ojj  ^  0,  022  “  ^^12  “  ®  ^  “  (/2/3  ^  “  (*^  2 

®22  ^  ‘^ll  “  ^12  “  ®  “  (/2/3)de^22’  ^  ~ 

(2.22) 

The  case  for  1-D  stress  at  an  arbitrary  orientation  is  dealt  with  below. 
Uniaxial  Specimens 

For  a  uniaxial  specimen  cut  from  an  orthotropic  material  at  an 
arbitrary  orientation  to  the  principal  axes,  the  state  of  stress  is 
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transformed  back  to  the  material  directions  as  follows 


^11  ~  OgCOS  0 
2 

^22  ® 

0^2  “  OgSinGcosS  (2.23) 


(See  Figure  13  for  an  illustration  of  the  typical  uniaxial  specimen.) 
The  incremental  strains  in  the  loading  direction  and  in  the  transverse 
direction  are 


dCg  =  de^jj^cos^e  +  de^22®^’'^®  2de^ ^2®^’^®^°®® 

4  2  2 

=  [rj^j^cos  6  +  2rj^2COs  0sin  0 

4  2  2 

+  sin  0  +  2r  sin  0cos  0]dX0 
66  0 

de^^Q  =  de^^j^sin^0  +  de^22‘^°®^®  "  2de^ ^2*^°®®®^’^® 

2  2 

=•  [r, ,  +  1  -  2r^,)sin  0cos  0 
11  66 

4  4 

+  r,„(sin  0  +  cos  0)ldXa„ 

IZ  u 


(2.24) 


The  effective  stress  and  the  incremental  effective  plastic  strain  are 
found  to  be 

—  4  4  2  2  ' 

a  =  /  y  [tjjCos  0  +  sin  0  +  2(r^^  +  rj^2)®i"  0cos  0]  Og 

d^  =  de^-// [r,  cos^0  +  sin^0  +  2{r  +  r.  )sin^0cos^6]'  (2.25) 

u  L  i  1  DO  \-Z 

The  ratio  do/de*^  can  be  derived  from  the  expression  for  dX.  Substituting 

i->n 

the  expressions  for  a  and  de  .  to  obtain 

u 
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dX  “  2  (d^/da)  (da/a) 


4 

r^jCos  6 


+  sin^0 


+  2(r 


66 


2  2 

+  rj2)sin  6cos  6 


de^e 

da„ 


dg^ 

^6 


(2.26) 


Hence  obtain 

da/d^  “  -2  [rj^jcos  0  +  sin  0 

+  2(rgg  +  rj2)sln^6cos^0]  (dag/de^^)  (2.27) 

The  plastic  Poisson's  ratio  effect  Is  defined  as  v^0»  takes  the  form 

-vP  =  =  ^^11  -*•  ^  ~  -»• 

®  rjj^cos^0  +  sln^O  +  2(rj2  +  rgg)sln^0cos^0  (2.28) 


Note  that  for  various  orientations  the  sine  or  cosine  terms  are  zero  and 
the  expression  for  v^q  simplifies,  l.e. 


0-0'* 

0  »  90“ 
0  =  45“ 


'12/^1 
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^11  '  ^^66  ^^12 
+  2^66  +  2rj2 


(2.29) 


These  relations  are  useful  because  they  give  a  general  Idea  of  how  the 
vP-  versus  0  plot  will  look  depending  on  the  r-values. 

EXPERIMENTAL  TESTING  AND  ANALYSIS 

By  testing  uniaxial  specimens  cut  at  various  angle  to  the  fibers, 
two  relevant  sets  of  data  are  taken  versus  the  applied  stress,  l.e., 
longitudinal  strain  and  transverse  strain.  The  total  strains  must  then 
be  separated  into  elastic  and  plastic  components. 
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Elastic  Characterization 


The  s-values  are  solved  for  by  curve-fitting  a  plot  of  the  elastic 
Poisson's  ratios  versus  the  fiber  orientation.  Note  that  in  order  to 
determine  the  three  s-values,  at  least  three  data  points  are  required. 
These  values  do  not  serve  to  characterize  the  material  alone;  as  was 
stated  before,  E2  =  da/de  =  E,  and  with  the  s-values  already  calculated 
the  relationship  exists  to  reduce  the  collected  E-  values  into  the 
effective  elastic  modulus,  E.  Therefore,  at  least  one  value  of  E^  is 
required  to  calculate  E.  Once  E2  is  obtained,  the  other  principal 
direction  elastic  constants  can  be  obtained  from  the  s-values  using  the 
relations  given  in  equation  (2.3). 

Note  that  it  is  perhaps  just  as  effective  to  determine  Ej^,  E2, 

Gj^2»  experimental  testing  as  described  in  Reference  [1]. 

However,  such  tests  presuppose  that  the  material  is  orthotropic,  while 
reducing  experimental  data  with  the  analysis  described  in  this  chapter 
effectively  checks  for  material  orthotropy. 

Plastic  Characterization 

The  plastic  data  collected  is  organized  into  two  forms.  We  have; 

Observe  that  E^  is  essentially  similar  to  an  elastic  modulus,  while 

0 

is  the  plastic  Poisson's  ratio.  It  should  be  noted  that  both  E^g 
and  may  be  nonlinear  functions  of  strain. 

Once  the  r-values  have  been  obtained,  it  is  then  possible  to  reduce 
the  stress  versus  plastic  strain  data  into  effective  stress  versus 
effective  plastic  strain.  If  the  r-values  are  good,  the  results  of  this 
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reduction  will  tend  to  collapse  the  stress-strain  data  together  so  as 
to  form  a  single  curve.  This  effective  stress-strain  curve  is  then 
curve-fitted  as  appropriate.  The  curve-fit  should  be  one  which  allows 
for  a  derivative  to  be  easily  taken.  It  is  this  derivative  (da/d?*) 
which  is  necessary  to  complete  the  analysis.  The  final  form  for  the 
orthotropic  plastic  characterization  is  given  below 


where 


’^ll  ^^12 


^^12  1 


0 

0 


0  2r 


66 


11 


22 


'12 


d?/da  is  found  from  the  curve-fit 


(de^/da) (da/a) 


(2.30) 


a  = 


3  /  2  2  2 

2  *^*‘ll‘^  +  a  +  2r,.,a,  ,a«.,  +  2r^^a 


22 


12“ir22 


66'^  12 


3 


do  -  2  a/0)((rj20,2  +  r22022)do,j 


+  (rjjOjj  +  022)d022  +  2-^66 '’l2'*‘’l2> 


The  values  of  t’j2»  are  found  from  the  curve-fit  of  the 

versus  9  data.  Note  that  again  three  values  of  are  required  in 

y 

order  to  determine  the  three  r-values,  and  at  least  one  set  of  a  versus 
e^  data  is  required  to  find  d?/da  . 


Limiting  Values 

Reference  [3]  gives  several  expressions  which  describe  limiting 
values  of  the  elastic  orthotropic  material  constants.  The  expressions 
are  derived  from  work  considerations  and  the  symmetry  of  the  compliance 
matrix.  Such  restrictions  are  used  to  check  experimental  data  to  see 
if  they  are  consistent  within  the  framework  of  the  elasticity  model. 
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Similar  limiting  cases  can  be  derived  for  the  r-values  determined  in 
the  plastic  analysis. 

From  equation  (2.2  4),  the  incremental  relation  for  xiniaxial  stress  is 

D  4  4  2  2 

de*^g  =  [rj^j^cos  9  +  sin  0  +  6cos  0]agdX 

Since  de^.  must  always  be  positive  for  a  positive  a.,  then 

y  y 

4  4  2  2  > 

r^j^cos  0  +  sin  0  +  "*■  0cos  0-0 

or  r^^j  =  0  and  rj2  +  =  0  (2.31) 

Now  we  let  r^j^  =0  to  obtain  the  following 

T\  9  9 

de^^g  =  [(1^21  ^  9cos  9 

4  4 

+  ®  0)l<^0dX 

=  [(1  -  2r ^0)sin^0cos^0 

+  r,_(sin^9  +  cos^0)]aQdX  (2.23) 

iz  y 


In  order  to  avoid  having  a  width  expansion  under  a  tensile  load,  then 


and 

hence 


or 


In  summary,  we  have  these  limiting  conditions  which  the  experimental 
data  should  obey 


and 


(2.33) 
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For  an  isotropic  material  it  is  evident  from  the  definitions  of  the 


r-values  that 


'll  ■  1 


'12  '  --5 


'66  - 


Therefore  the  r-values  for  the  isotropic  case  are  within  the  limitations 
of  the  mathematical  modeling  of  the  plasticity.  In  Chapter  3  the  r-values 
obtained  from  the  analysis  on  the  experimental  data  will  be  checked 
against  these  limiting  conditions  as  well. 


COMPUTER  MODELING 

The  equations  for  the  elastoplastlc  characterization  were  put  into 
a  program  to  generate  the  stress-strain  curves  of  a  uniaxial  specimen 
at  an  arbitrary  orientation.  The  program  reads  in  the  material  properties 
and  fiber  orientation,  calculates  material  stresses  by  an  angle  trans¬ 
formation,  then  calculates  elastic  strains.  A  yield  criterion  based  on 
the  effective  stress  is  used,  and  if  the  specimen  has  yielded  the  work 
hardening  model  is  employed  to  calculate  the  plastic  strains.  The  sum 
of  the  elastic  strain  and  plastic  strain  is  output  versus  the  uniaxial 
stress,  for  both  longitudinal  and  transverse  strain.  The  work  hardening 
model  used  is  the  expression  chosen  to  curve-fit  the  effective  stress- 
strain  curve  of  the  experimental  data.  The  program  is  called  UNIAX.F 


Parametric  Study 

Example  rxms  were  conducted  to  gain  an  understanding  of  orthotropic 
plasticity.  The  work  hardening  curve  used  was  a  simple  three-parameter 
model  as  given  below: 


-  (Oy/A)® 


(2.34) 
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The  values  used  for  the  parameters  were  fixed  at  A  =  50000  psi,  B  =  5, 
and  cr  =  20000  psi.  Note  that  it  is  this  value  of  a  which  is  also  used 

y  y 

in  the  yield  criterion. 

The  study  focused  on  two  areas;  isotropic  versus  anisotropic  and 
elastic  versus  plastic.  Three  orientations  were  used  in  each  case: 

0®,  45®,  and  90®.  Obviously,  for  the  isotropic  elastic  and  isotropic 
plastic  case,  the  stress-strain  history  is  expected  to  be  invariant 
with  respect  to  specimen  orientation.  The  program  proved  capable  of 
generating  this  trivial  case. 

The  nex"  case  was  to  check  for  anisotropic  elastic  and  isotropic 
plastic  behavior.  To  accomplish  this,  the  following  values  were  used: 

Ej^  =  20  msi  =  2  msi  ~  ~ 

’^ll  "  ^  ’^12“  ^66  “ 

The  plots  of  stress  versus  longitudinal  strain  and  transverse  strain 
versus  longitudinal  strain  are  shown  in  Figures  23  and  24  respectively. 
The  points  to  observe  in  Figure  23  are  that  the  elastic  slopes  are  indeed 
dependent  on  the  fiber  angle,  and  the  yield  is  independent  of  the  fiber 
angle.  It  is  not  obvious  from  this  plot  that  the  work  hardening  curves 
obey  the  same  law  for  each  orientation,  but  it  is  true  in  this  case. 

Note  that  the  material  yields  at  a  value  less  than  =  20000  psi.  This 
is  due  to  the  definition  of  the  effective  stress  which  is  used  in  the 


yield  criterion.  The  true  value  of  yield  is 


for  an  isotropic  plastic  material. 


/To  =  16330  psi 
3  ^ 


In  contrast  to  the  previous  case,  the  parameters  for  an  isotropic 


elastic  and  anisotropic  plastic  material  were  input.  (Again,  r-values 
were  chosen  which  were  within  the  limiting  conditions  for  an  orthotropic 


material.)  These  values  were 


Ej^  =  20  msi  E2  =  20  msi  “  7.6923  msi  Vj2  “  *3 

■^11  '  ^2  ’  '66  ■ 

The  plots  for  this  case  arc  shown  in  Figures  16  and  17.  Note  that  the 
elastic  curves  are  all  independent  of  the  load  angle  except  that  for 
each  orientation  the  yield  point  is  different.  Indeed,  the  0®  specimen 
has  not  yielded  on  this  plot,  "^is  variation  of  yield  with  fiber  direction 
is  a  product  of  the  anisotropic  r-values.  Although  it  is  not  obvious 
from  Figure  16,  the  work  hardening  curves  are  also  different  for  each 
orientation. 

Combining  the  above  two  cases  into  an  anisotropic  elastic  and 
anisotropic  plastic  case  Incorporates  the  features  of  elastic  slopes, 
yield  points,  and  work  hardening  curves  which  vary  with  fiber  orientation. 
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Chapter  3  MATERIAL  CHARACTERIZATION 


The  following  experimental  tests  for  the  characterization  of  the 
elastic  and  plastic  orthotropy  were  designed  according  to  the  requirements 
set  forth  in  Chapter  2. 

The  tests  were  intended  primarily  to  obtain  the  work  hardening 
curves  i.e.  designed  for  large  strains,  but  elastic  properties  obtained 
were  of  sufficient  quality  to  perform  the  elastic  characterization. 

TEST  SPECIMENS  AND  DATA  COLLECTION 

The  tensile  specimens  for  the  testing  were  cut  from  a  second  boron/ 
aluminum  panel  of  Identical  manufacture  to  the  first.  The  cutting  scheme 
for  specimens  17  through  24  appears  in  Figure  18.  Figure  19  shows  the 
typical  test  specimen.  Each  specimen  was  nominally  nine  inches  long 
with  aluminum  tabs  applied  with  cyanoacrylate  glue,  leaving  a  test 
section  of  nominally  six  inches.  The  specimen  widths  were  all  approxi¬ 
mately  one-half  inch.  This  gave  a  length-to-width  ratio  of  12  which  was 
considered  to  be  adequate  for  the  fiber  orientations  chosen.  The  exact 
dimensions  and  the  fiber  orientations  for  the  specimens  are  summarized 
in  Table  3.  Loading  was  uniaxial  tension,  with  fixed-grip  end  conditions. 

To  measure  the  longitudinal  and  transverse  strains  the  clip  gage  described 
in  Chapter  1  was  used.  The  output  from  this  gage  was  digitized  and 
stored  on  floppy  disk.  The  schematic  for  the  test  setup  and  the  data 
acquisition  system  is  shown  in  Figure  20.  The  MTS  was  operated  under 
head-displacement  control  to  observe  the  extended  yielding  of  the  spe-iimens. 
Loading  was  controlled  manually,  and  after  each  increment  the  Daytronics 
was  cycled  to  record  the  strains  and  load  simultaneously.  The  data  was 
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stored  on  the  Apple  III  computer  by  means  of  a  small  program  which  read 
from  the  RS232  port  of  the  Daytronics  and  wrote  directly  to  disk. 

The  Clip  Gage 

Both  longitudinal  arms  of  the  clip  gage  were  wired  into  one  channel 
to  complete  a  full-bridge,  and  the  transverse  strain  arm  was  wired  in 
as  a  half -bridge.  A  variable  potentiometer  circuit  was  constructed 
to  step  down  the  load  voltage  from  the  MTS  machine.  Care  was  taken  not 
to  have  a  mismatch  of  Impedance  when  designing  the  potentiometer.  The 
resulting  circuit  was  capable  of  reducing  the  output  voltage  by  a  factor 
of  one  thousand  so  that  the  Daytronics  could  accept  the  load  voltage 
into  one  of  the  strain  conditioners. 

Since  the  clip  gage  was  designed  for  one-inch  wide  specimens,  a 
recalibration  of  it  for  the  narrow  specimens  was  conducted.  Figures  21 
and  22  show  the  comparison  nm  on  a  half-inch  wide  piece  of  aluminum 
with  well  known  material  properties.  Two  strain  gages  (front  and  back) 
were  mounted  on  the  specimen  to  provide  a  standard  to  which  the  clip  gage 
was  compared.  The  deviations  of  the  data  from  the  45®  line  as  shown 
in  Figure  22  was  due  to  the  extended  pins  used  on  the  clip  gage  to 
secure  it  to  the  narrow  specimen.  These  longer  pins  caused  an  increased 
moment  arm  which  twisted  the  gage  arm,  effectively  decreasing  the  observed 
strain.  A  calibration  factor  of  1.5504  was  used  for  all  subsequent  tests 
on  half -inch  wide  specimens. 

Data  Reduction 

The  results  of  the  tensile  tests  are  shown  grouped  together  in 
Figures  23  through  26.  Each  plot  combines  the  data  for  all  orientations. 
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Because  of  the  different  orders  of  magnitude  of  ultimate  stresses  and 
ultimate  strains  between  them,  little  detail  can  be  seen  in  Figures  23 
and  25.  For  this  reason  Figures  24  and  26  give  an  enlarged  view  of  the 
yield  curves.  As  can  be  seen  from  these  plots,  the  tests  give  similar 
results  for  the  0®,  10®,  and  90®  orientations  as  compared  with  the  results 
presented  in  Chapter  1.  More  accurate  Information  concerning  the  trans¬ 
verse  strain  was  recorded  in  these  tests,  however.  As  seen  best  in 
Figure  26  the  elastic  and  plastic  Poisson's  ratios  range  from  very  low 
in  the  90®  specimen  to  very  high  in  the  10®  and  20®  specimens.  The 
elastic  and  plastic  constants  calculated  from  these  tests  are  tabulated 
for  each  orientation  in  Table  4. 

For  each  orientation  the  first  step  was  to  obtain  the  elastic 
modulus.  With  this  value  the  elastic  strains  were  subtracted  from 
the  total  strains  to  give  the  plastic  strains  as 


(3.1) 


These  two  steps  were  accomplished  with  one  computer  program.  Input  into 
the  program  was  the  stress  versus  total  strain  file  along  with  the 
number  of  data  points  selected  to  represent  the  total  extent  of  the 
elastic  (or  linear)  portion  of  the  curve.  The  program  output  the  Young's 
modulus  and  the  elastic  Poisson's  ratio.  It  then  subtracted  the  elastic 
strains  according  to  the  above  equations  and  wrote  the  plastic  strain 
data  into  a  file.  A  second  program  was  then  used  to  obtain  the  plastic 
Poisson's  ratios  from  these  files  of  plastic  longitudinal  strain  versus 
plastic  transverse  strain. 
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ELASTIC  ORTHOTROPIC  AMALYSIS 


A  plot  of  elastic  Poisson's  ratio  versus  fiber  orientation  is  given 
in  Figure  27.  (Note  that  the  data  reported  from  Chapter  1  shows  consi¬ 
derable  scatter.  Also,  the  value  of  Poisson's  ratio  obtained  from  the 
65“  specimen  is  considered  to  be  in  error.)  This  data  was  curve-fit 
to  obtain  the  values  of  Sj^2  boron/altinimum. 

Figure  28  is  a  plot  of  the  elastic  moduli  versus  fiber  orientation. 
The  solid  line  is  a  least-squares  fit  of  the  data.  The  s-values  obtained 
from  Figure  27  were  used  to  reduce  both  the  data  points  and  the  curve-fit 
into  effective  elastic  moduli  values.  The  dashed  line  is  the  reduced 
curve-fit,  and  it  is  a  constant  with  fiber  direction  as  expected,  with 
a  value  of  19.1  msi.  Hence,  the  following  values  were  obtained  for  the 
principal  material  directions 

Ej  =  29.4  msi  E2  =  19.1  msi 

Vji2  =  *169  Gj^2  “  msi 

These  values  compare  favorably  with  those  reported  in  other  literature. 

PLASTIC  ORTHOTROPIC  ANALYSIS 

Figure  29  shows  the  plot  of  the  plastic  Poisson's  ratio  values 
obtained  for  the  various  orientations.  The  figure  also  shows  two  curve- 
fit  lines  for  the  data.  The  dashed  line  is  a  least-squares  curve  fit 

which,  incidently,  intercepts  the  y-axis  at  a  value  of  =  4.5. 

u 

The  r-values  obtained  from  this  curve-fit  were  used  to  calculate  effective 
stress  and  effective  strain  plots  from  the  plastic  data  and  are  plotted 
together  as  shown  in  Figure  30.  The  effect  of  missing  the  0“  data  with 
the  curve-fit  can  be  seen  quite  clearly.  All  other  data  sets  fall 
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fairly  close  together.  As  was  stated  In  Chapter  2,  obtaining  this  master 
curve  was  the  primary  goal  of  the  analysis  so  that  it  may  be  curve-fit 
to  model  the  material  behavior.  Because  of  the  large  error  bounds  on 
the  plastic  Poisson’s  ratio  data  (perhaps  as  much  as  ±.l)  it  was 
questionable  whether  the  r-values  given  by  the  least-squares  method  were 
optimal  for  obtaining  a  master  curve  of  effective  stress  and  strain.  A 
small  program  was  written  which  allowed  the  r-values  to  be  input  manually 
to  observe  how  the  resulting  output  modeled  the  data.  Working  by  trial 
and  error  the  solid  line  curve  shown  in  Figure  29  was  obtained  for  the 
following  r-values 


r,i  -  .017 

'12  °  --“I 


■^66  “  I-® 


(3.2) 

Observe  that  these  values  are  within  the  limiting  conditions  stated  in 
Chapter  2.  Perhaps  these  values  are  not  optimal  either,  but  the  resulting 
master  curve  is  superior  to  the  one  obtained  from  the  least-square 
r-values.  For  the  r-values  given  in  (3.2)  the  master  curve  in  Figure  31 
was  obtained.  This  plot  also  shows  the  parametric  curve  which  was  used 
to  fit  the  data.  It  was  found  that  the  functional  form  from  (2.34)  was 
sufficient  for  modeling  the  data, 

0  ; 


eP  = 


—  <  — 
a  =  a 


(a/A)®  -  (Oy/A)® 


A  =  64000  psi 


B  =  5.8 


—  >  — 
a  =«  o 


a  =  14500  psi 

y 


(3.3) 


This  is  the  functional  form  of  the  work  hardening  model  which  is  incor¬ 
porated  into  the  computer  programs,  along  with  its  derivative  which  is 
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easily  obtained.  Using  this  master  curve,  the  stress  versus  plastic 
strain  curves  for  any  orientation  can  be  generated. 

Coiid)lning  the  elastic  and  plastic  strains,  the  complete  stress- 
strain  histories  for  any  fiber  orientation  can  be  produced.  Figures  23 
through  26  are  plots  of  the  generated  stress-strain  curves  superliiq)osed 
on  the  actual  experimental  data.  It  can  be  observed  that  the  generated 
plots  capture  the  essence  of  the  material  behavior. 

DISCUSSION 

It  is  obvious  from  the  figures  of  the  stress-strain  data  collected 
that  some  specimens  displayed  considerable  plastic  strain  while  others 
failed  at  relatively  low  values  of  plastic  strain.  This  gives  rise  to 
the  question  of  damaged  specimens.  While  damage  due  to  cutting  is  of 
concern  especially  as  the  specimens  become  thinner,  this  is  only  a 
secondary  consideration  compared  to  the  damage  which  may  be  caused  by 
gripping  the  tab  ends  with  too  much  pressure.  Specimen  16  for  exaii;>le, 
which  was  the  first  half -inch  wide  specimen  tested,  failed  at  a  stress 
of  141  ksl,  not  nearly  the  200  ksi  as  expected.  When  the  specimen  was 
removed  from  the  MTS  grips  it  was  observed  that  the  material  under  the 
tabs  had  been  deformed  by  excessive  grip  pressure.  Specimen  21  was  cut 
to  perform  the  test  over;  its  data  and  not  that  from  specimen  16  was 
used  in  the  analysis.  More  care  was  taken  in  subsequent  tests  to 
avoid  end  damage,  yet  several  of  the  specimens  fractured  near  the  tab 
ends  regardless.  A  summary  of  the  fracture  modes  and  the  ultimate 
strains  from  the  specimens  is  presented  in  Table  5.  Under  Failure 
Mode,  the  term  'end'  means  that  the  specimen  fractured  at  the  tab. 


35 


and  'middle*  indicates  that  the  fracture  occurred  away  from  the  tab. 

From  this  table,  it  is  obvious  that  there  is  a  correspondence  between 
low  values  of  longitudinal  failure  strains  and  the  fractures  which 
occurred  at  the  tab  ends. 

Concerning  the  analysis,  it  should  be  noted  that  the  material  at 
a  0®  orientation  behaves  in  a  singular  manner.  All  other  orientations 
yielded  stress-strain  curves  typical  of  most  ductile  metals.  Only  the 
0®  specimens  behaved  as  a  brittle  material.  Because  it  was  noted  that 
the  specimens  showed  no  residual  strain  even  when  stressed  nearly  to 
failure,  the  material  in  this  orientation  could  be  considered  linear 
elastic  up  to  failure.  In  this  case,  there  would  be  no  plastic  strain 
data  and  consequently  an  indeterminate  plastic  Poisson's  ratio  for  the  0® 
specimen.  The  modeling  parameters  as  used  here  will  predict  some 
plasticity  for  a  0®  specimen.  However,  this  is  thought  to  be  of  little 
consequence  since  in  a  multlaxlal  stress  state  the  other  components  of 
stress  will  dominate  the  yielding. 
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Chapter  4  STRESS  ANALYSIS  OF  CRACKS 


Determining  the  stress  and  strain  distribution  around  the  crack 
tip  is  of  primary  importance  In  this  chapter.  This  will  be  done 
experimentally  by  mounting  very  small  strain  gages  close  to  the  crack. 

An  approximate  analysis  will  be  done  based  on  an  elastic  analysis  of 
a  crack  in  an  infinite  orthotropic  sheet. 

EXPERIMENTAL  STUDY 

In  order  to  examine  the  plastic  zone  at  the  crack  tip  in  boron/ 
aluminum,  experiments  were  run  with  strain  gages  mounted  at  the  crack 
tip.  One  test  was  on  a  crack  along  the  fibers  perpendicular  to  the 
load,  and  the  other  for  a  crack  along  the  fibers  at  a  60°  angle  to  the 
load.  The  data  collected  was  reduced  to  the  form  of  applied  remote 
stress  versus  strain.  The  objective  of  this  series  of  tests  was  to 
examine  the  plastic  zone  sizes  and  to  observe  the  plastic  zone  growth 
during  a  1.5  overload  cycle. 

Tapered  Specimens 

A  single  case  of  plastic  zone  growth  was  designed  In  the  form  of  a 
tapered  specimen  (specimen  25)  to  test  the  setup  and  material  modeling. 
Figure  48  shows  the  dimensions  and  the  gage  locations  on  the  specimen. 
Quarter-inch  grid  strain  gages  were  used  for  this  test.  The  data 
acquisition  setup  is  the  same  as  shown  In  Figure  20. 

Because  of  the  special  geometry  required  of  the  specimen,  the 
diamond  saw  could  not  be  used.  The  specimen  was  cut  by  electric  discharge 
machining  (EDM) .  The  cut  edge  was  found  to  be  of  superior  quality  to 
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that  of  the  saw  cut  edge.  However,  although  the  specimen  is  immersed 
in  oil  during  cutting,  it  is  questionable  whether  or  not  the  material 
edges  are  subject  to  excessive  heating  during  the  electric  discharge 
machining. 

The  data  from  the  tapered  specimen  can  be  found  in  Figures  33  and  34, 
for  the  longitudinal  gages  only.  The  data  in  Figure  33  is  arranged  as 
applied  tensile  load  versus  strain.  As  expected,  the  gage  at  the  narrow 
end  (gage  1)  shows  signs  of  yielding  first,  and  the  gage  at  the  wider  end 
(gage  6)  displays  yielding  last.  Thus,  a  progression  of  plastic  zone 
is  observed.  (Gages  5  and  7  were  placed  front  and  back  in  the  same 
location  to  check  for  bending.  Both  gages  show  similar  strains.)  The 
data  was  also  reduced  to  the  form  of  net  section  stress  versus  strain  to 
see  if  the  data  generally  collapsed  onto  the  same  line.  Figure  34 
shows  the  plot  of  this  data,  as  well  as  results  from  a  finite  element 
modeling.  The  model  data  predicts  larger  strains  for  a  given  stress  level 
or  conversely,  the  tapered  specimen  appears  to  be  more  constrained  than 
either  the  uniaxial  specimen  or  the  finite  element  model.  The  reason 
for  the  discrepancy  appears  to  be  due  to  the  long  high  strength  fibers 
extending  from  the  elastic  region  into  the  plastic  region.  The  main 
effect  of  these  fibers  is  to  inhibit  the  shear  plastic  flow  and  this  is 
corrobarated  by  Figure  34  which  also  shows  a  finite  element  modeling  but 
with  the  shear  coefficient  changed  from  1.8  to  1.5,  and  the  agreement 
is  very  close. 

The  above  result  is  important  when  considering  the  modeling  of  the 
crack  tip  in  the  orthotropic  material  because  there  is  a  substantial 

elastic  constraint  on  the  strains  in  the  region  of  local  yielding. 
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In  other  words,  the  yielded  material  is  not  allowed  to  flow  as  freely 
due  to  the  strong  fibers  protruding  into  the  plastic  zone.  It  is  clear 
that  in  order  to  account  for  this  constraining  effect  either  a  finite 
element  program  which  models  the  actual  fibers  is  required,  or  some 
testing  scheme  for  obtaining  the  'r*  parameters  under  constrained 
conditions  be  devised.  The  former  seems  too  cumbersome,  but  the  latter 
has  potential  if  a  tapered  specimen  is  used  in  the  fundamental  test. 

Strain  Gage  Strips 

The  choice  of  resistance  gages  for  measuring  the  strains  at  the 
crack  tip  was  based  on  the  need  for  ease  and  reliability  of  data  acqui¬ 
sition.  It  was  understood  that  the  strain  gages  would  be  operating  under 
the  worst  of  circumstances:  A  steep  strain  gradient  and  the  size  of  the 
plastic  zone  required  very  small  gages,  whereas  normally  the  inhomogeneity 
would  require  large  gages  to  average  out  the  strains  over  a  greater 
area.  The  gages  purchased  from  Mlcro-Weasurements ,  Inc.  were  in  the  form 
of  a  strip  of  ten  gages,  of  alternating  grid  directions,  type  EA-13-031MH-120. 
The  grid  sizes  of  the  gages  were  each  .031"  across,  and  the  entire  strip 
was  on  the  order  of  only  .7"  long.  The  strip  was  cut  into  two  sections, 
one  of  six  gages  and  one  of  four  gages.  For  each  of  the  two  tests  these 
strips  were  placed  on  the  X  and  Y  axes  relative  to  the  crack  tip.  The 
sequence  for  mounting  the  gages  was  to  first  initiate  a  crack,  then  the 
gages  were  moxmted  a  short  distance  away  from  the  crack  tip.  The  specimen 
was  fatigued  under  a  constant  amplitude  cyclic  stress  until  the  tip  was 
close  to  the  first  gage  grid.  (Note  that  it  aided  the  positioning  of 
the  gages  to  know  that  the  crack  would  consistently  travel  in  the  fiber 
direction. ) 
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Study  of  Transverse  Crack 


Specimen  26  was  prepared  such  that  the  applied  stress  would  be 

perpendicular  to  the  crack.  The  strain  gages  were  moimted  as  shown  in 

Figure  35  and  the  crack  was  grown  into  position  so  that  during  the  over¬ 
load  cycle  the  gages  were  at  the  (X,Y)  coordinates  given. 

The  plots  for  gages  1,  3,  5,  and  the  remote  gage  are  shown  together 
in  Figure  36.  (Note  that  gage  1,  being  right  next  to  the  crack  tip, 
failed  during  the  test.)  The  amount  of  plasticity  recorded  corresponds 
to  the  distance  from  the  crack  tip.  It  is  expected  that  the  yield  stress 
at  each  gage  should  also  correspond  to  the  distance  from  the  crack  tip. 
Since  it  is  evident  that  the  apparent  yield  stresses  for  gages  1,  3,  and 
5  are  all  at  9400  psi  (which  is  the  maximum  cyclic  stress  prior  to  the 
overload) ,  it  can  be  concluded  that  the  material  at  these  locations 

yielded  prior  to  the  1.5  overload.  Observe  also  that  the  remote  gage 

shows  some  residual  strains  following  the  overload  cycle.  This  indicates 
that  the  entire  specimen  has  yielded  in  post  overload  to  a  small  degree. 

Study  of  Oblique  Crack 

Specimen  27  was  cut  so  that  the  crack  would  be  at  an  angle  of  60“ 
to  the  applied  stress.  Figure  37  is  a  schematic  of  the  specimen,  giving 
geometry  and  (X,Y)  gage  coordinates  prior  to  the  1.5  overload  application. 

Figure  38  shows  the  data  for  gages  1,  3,  5,  and  the  remote  gage. 

The  maximum  stress  of  the  constant  amplitude  cycling  was  made  lower  for 
this  test  to  avoid  yielding  the  entire  specimen  as  in  the  previous  test. 

It  can  be  observed  from  the  remote  gage  that  there  is  very  little  residual 
strain  after  the  overload.  Again,  as  expected,  the  amount  of  residual 
strains  in  gages  1,  3,  and  5  correspond  to  the  distance  from  the  crack  tip. 
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This  time,  however,  the  apparent  yielding  occurs  at  higher  stresses  the 
further  the  gage  is  from  the  crack  tip.  This  indicates  that  the  gages 
have  recorded  an  expansion  of  the  plastic  zone  during  the  overload. 

APPROXIMATE  ELASTIC-PLASTIC  ANALYSIS  OF  CRACKS 

A  computer  program  was  developed  in  order  to  obtain  an  estimate  of 
the  stresses  around  the  crack  tip  in  B/Al.  This  program  was  used  to 
predict  plastic  zone  shapes  and  sizes,  and  to  give  stress  and  strain 
predictions  at  arbitrary  points  around  the  crack  tip. 

A  Crack  in  an  Infinite  Orthotropic  Sheet 

For  an  approximte  crack  analysis,  a  program  was  written  which 
made  use  of  the  solutions  for  the  elastic  stresses  (mode  I  and  mode  II) 
in  an  orthotropic  sheet.  Incorporated  in  the  program  is  the  effective 
stress  yield  criterion  and  the  work  hardening  model.  The  expressions 
for  the  stresses  were  derived  from  the  elastic  solution  of  the  stress 
distribution  around  an  elliptical  hole  [Reference  5].  The  tractions  at 
the  surface  of  the  ellipse  are  required  to  be  zero  and  the  stresses  at 
infinity  are  required  to  be  the  applied  stress.  The  total  stresses  are 
the  sums  of  two  solutions,  as  shown. 

2 

a  =  a  '  +  a  cos  'V 

XX  «> 

2 

a  =  a  '  +  a  sin  ¥ 

y  y  00 

a  =  a  '  +  a  sinYcosY  (A.l) 

xy  xy  <» 

where  the  sign  convention  and  notations  are  shown  in  Figure  39* 

The  values  o  ' ,  a  ’ ,  and  a  '  are  stresses  resulting  from  the  tractions 
X  y  ’  xy  ® 

on  the  ellipse.  The  crack  solution  is  obtained  by  letting  the  minor  axis 
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of  the  ellipse  be  equal  to  zero.  In  this  limit,  the  stresses  are 
related  to  the  stress  functions  and  4>2  by  the  equations 


where 


and 


and 


a  '  =  2Re{$,  + 

y  12 

-  a^(l/pj) 

AGO 

*  —  {asin'J' (cosf  +  ^ 

^21  ~  ~  {asinf  (cosH'  +  Uj^sin'l')  } 


The  singular  stresses  could  be  derived  from  these  expressions  by  taking 
the  limit  of  the  distance  from  the  crack  tip  to  be  much  less  than  the 
crack  length  (r  <<  a).  In  this  case,  -»■  1  and 


Zj^  =  X  +  y^y  »  a  +  rcos9  +  y^sinG 


so  that 


-  a^  v4r(cos9  +  y^sin0) 


This  gives  the  near-tip  solution  with  ^  singular  stresses  at  the  crack 
tip  for  mode  I  and  mode  II.  The  program  ,  however,  does  not  make  this 
approximation,  and  will  be  referred  to  as  ELLPLAS.F. 

Using  the  r-values  and  the  elastic  constants  obtained  for  boron/ 
aluminum,  the  plastic  zone  shapes  were  plotted  for  fibers  running 
parallel  to  the  crack  and  perpendicular  to  the  crack.  These  plots  are 
found  in  Figures  40  and  41.  By  comparison,  the  program  was  also  made  to 
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output  the  plastic  zone  shape  for  an  isotropic  material  as  shown  in 
Figure  42.  Each  zone  was  calculated  for  the  same  remote  stress  (applied 
perpendicular  to  the  crack)  and  the  same  crack  length.  The  zone  shape 
resulting  from  the  fibers  being  perpendicular  to  the  crack  is  worth  noting. 

The  effect  of  changing  the  loading  angle  with  respect  to  the  crack 
can  be  observed  in  Figures  43  through  45.  In  each  plot  the  fibers  are 
parallel  to  the  crack.  Starting  with  the  load  perpendicular  to  the 
crack,  it  is  shifted  to  60®  and  then  to  30°.  The  zone  shape  changes  as 
mode  II  becomes  predominant,  and  the  size  of  the  zone  shrinks  as  the 
load  approaches  0®  (parallel  to  the  crack)  where  no  singular  stresses 
would  be  observed. 

Reconsidering  the  Strain  Gage  Results 

An  estimate  of  the  plastic  zone  corresponding  to  the  transverse  crack 
test  can  be  seen  in  Figure  41.  The  correct  crack  length  of  .31"  was 
used,  and  of  the  two  remote  load  levels,  the  9400  psi  zone  size  corresponds 
to  the  maximum  cyclic  stress.  The  load  level  of  14100  psi  corresponding 
to  the  peak  overload  stress  was  not  used  to  give  a  zone  size  because  it 
was  found  to  yield  the  entire  specimen.  To  better  illustrate  the 
problem,  the  same  program  was  used  to  plot  the  effective  stress  versus 
the  distance  from  the  crack  tip.  Figures  46  and  47  are  for  the  X  and  Y 
axes  respectively,  and  each  plot  also  shows  the  locations  of  the  strain 
gages.  Where  the  effective  stress  line  crosses  the  effective  yield 
stress  line  is  the  boundary  for  the  plastic  zone.  From  the  figures  it 
can  be  observed  that  only  half  of  the  strain  gages  are  predicted  to  be 
in  the  plastic  zone  for  the  stress  level  prior  to  the  overload.  This 
Indicates  that  the  approximation  underestimates  the  plastic  zone  size. 


43 


since  the  experimental  results  show  that  the  material  at  each  gage  location 
has  yielded  before  the  overload  cycle.  The  prediction  of  a  smaller  zone 
size  is  to  be  expected,  since  the  approximation  only  gives  an  elastic 
solution;  a  redistribution  of  the  stresses  due  to  the  yielded  material 
would  result  in  larger  predicted  zone  sizes.  This  effect  is  enhanced  further 
since  the  real  specimen  has  a  finite  width.  For  the  peak  overload  stress 
the  approximation  predicts  that  the  whole  specimen  yields,  as  can  be  seen 
in  the  two  figures. 

An  estimate  of  the  plastic  zones  for  the  oblique  crack  test  is 
shown  in  Figure  44.  Again,  to  better  illustrate  the  effective  stress 
distribution.  Figures  48  and  49  were  produced  similar  to  Figures  46  and  47. 
Unlike  the  previous  test,  the  plastic  zone  has  a  finite  size  for  the  peak 
overload  stress,  yet  the  estimates  show  that  each  gage  should  have  detected 
some  yielding  during  the  overload  cycle. 

DISCUSSION  OF  PLASTICITY  EFFECTS  AND  FATIGUE 

The  major  point  observed  in  the  fatigue  tests  mentioned  in  Chapter  1 
are  the  tendency  for  the  crack  to  always  grow  along  the  fibers.  The  tests 
conducted  in  this  chapter  offer  some  insight  into  these  phenomena. 

Regarding  the  crack  growth  direction,  the  plastic  zone  shapes 
predicted  by  the  approximate  analysis  give  an  indication  of  why  it  was 
found  to  be  difficult  to  grow  fatigue  cracks  through  the  fibers  in  B/Al. 
Figure  40,  which  shows  the  plastic  zone  extending  primarily  perpendicular 
to  the  crack  tip,  indicates  that  the  crack  would  be  prone  to  turning  and 
travelling  along  the  fibers. 
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CONCLUSIONS 


The  clip  gage  used  in  this  work  was  found  to  be  a  good  way  to 
determine  the  longitudinal  and  transverse  strains  in  a  coupon  specimen, 
but  it  is  recommended  that  specimen  dimensions  be  confined  to  those 
for  which  the  gage  is  designed  to  avoid  recalibration  inaccuracies. 

For  cutting  the  specimens  the  diamond  saw  was  adequate  as  long  as 
the  cuts  were  straight  and  the  specimens  were  not  required  to  have  sharp 
notch  cuts.  The  electric  discharge  machining  (EDM)  technique  is  capable 
of  cutting  any  shape  of  specimen  providing  an  electrode  of  that  shape 
can  be  made.  The  cutting  is  relatively  fast  and  flawless,  although  not 
necessarily  as  accurate  with  regard  to  specimen  fiber  orientation. 

Off-axis  specimens  should  be  cut  with  sufficient  length-to-width 
ratios  to  avoid  end  effects  if  uniaxial  testing  is  intended.  As  the  angle 
between  the  loading  direction  and  the  fiber  direction  becomes  smaller, 
there  is  a  need  for  larger  length-to-wldth  ratios.  A  ratio  of  12  was 
found  to  be  adequate  up  to  an  angle  of  20".  At  10“  the  end  effects  begin 
to  be  noticed,  as  indicated  by  the  data  presented  in  Figure  28.  The 
elastic  moduli  measured  for  the  10“  specimens  (both  specimen  5  and  specimen 
22)  were  consistently  higher  than  the  moduli  calculated  for  the  0“ 
specimens.  In  Reference  [1]  it  is  stated  that  a  length-to-wldth  ratio 
of  14  or  higher  be  used  at  a  fiber  orientation  of  10“ . 

From  the  discussion  in  Chapter  1  it  was  found  that  the  10“  off-axis 
tensile  specimen  is  not  optimal  for  determining  the  Intralaminar  shear 
of  a  mildly  orthotropic  material  such  as  boron/aluminum. 

The  elastic  orthotropic  characterization  given  in  Chapter  2  requires 
several  tests  (4  at  least) .  It  is  recommended  that  a  specimen  of 
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suitable  orientation  (determined  from  the  normalized  strain  versus  fiber 
angle  curve)  be  used  if  an  off-axis  test  is  to  be  performed. 

Without  doubt,  the  uniaxial  behavior  of  boron /aluminum  can  be  described 
within  an  anisotropic  elastic  plastic  formulation.  Two  crucial  questions 
remain  to  be  answered.  First,  what  is  the  behavior  of  the  material 
under  complex  loading  (Including  coii5)resslve  stresses.)  Secondly, 

(and  perhaps  more  importantly)  to  what  extent  does  the  inhomogeneity 
of  the  material  affect  the  plastic  flow  in  situations  of  constrained 
yielding?  Both  of  these  questions  can  be  clarified  only  after  more 
complicated  testing.  And  to  design  such  tests  will  require  a  finite 
element  program  capable  of  solving  orthotropic  elastic  plastic  problems. 

The  application  of  strain  gage  strips  to  determining  the  plastic 
strain  behavior  in  the  vicinity  of  cracks  appears  very  promising. 
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Table  1.  A  con^jrehensive  listing  of  B/Al  specimens  cut  from  panel 
1  as  shown  in  Figure  1. 


Specimen  and 
Fiber  Angle 

Dimensions  in  inches 
(t  X  W  X  1) 

Notch 

Cut,  (in) 

Use 

1 

0® 

.077  X  3.0  X  9.5 

.166 

.315 

.74 

Fatigue 

la 

0“ 

.077  X  1.92  X  9.5 

.74 

Fatigue 

lb 

0“ 

.077  X  1.05  X  9.5 

.122 

.281 

Fatigue 

2 

15“ 

.077  X  1.65  X  10. 

.465 

Fatigue 

3 

0“ 

.077  X  1.517  X  9.5 

.612 

Fatigue  + 

1.5  Overload 

4 

0“ 

.077  X  1.51  X  9.5 

.598 

Fatigue  + 

1.5  Overload 

5 

o 

O 

.077  X  1.005  X  8.6 

— 

6 

0“ 

.077  X  1.00  X  6.5 

— 

7 

0“ 

.077  X  1.011  X  10. 

— 

Material 

8 

90“ 

.0775  X  1.00  X  9.8 

— 

Properties 

9 

90“ 

.0775  X  1.006  X  9.8 

— 

10 

VO 

o 

o 

.077  X  3.08  X  9.5 

.148 

Fatigue 

11 

90“ 

.077  X  3.08  X  9.5 

.150 

Fatigue  + 

1.5  Overload 

12 

O 

O 

Ov 

.077  X  3.061  X  9.7 

.136 

Fatigue  + 

1.5  Overload 

13 

o 

O 

o^ 

.077  X  3.07  X  9.7 

.142 

Fatigue  + 
Fracture 

14 

0“ 

.078  X  1.52  X  9.5 

.290 

Fracture 

15 

0“ 

.077  X  1.51  X  9.5 

.571 

Fatigue 

Threshold 
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Table  2.  Mechanical  properties  of  B/Al  (tensile,  room  temperature) 


EXPERIMENTAL 

Elastic  Constants 


^1 

30.3  msl,  32.6  msl 

CM 

20.3  msl,  20.8  msl 

^12 

7.55  msl 

''l2 

.29 

^^21 

.21,  .37 

Ultimate 

Stresses 

204  ksl,  222  ksl 

^^2 

26.5  ksl,  27.0  ksl 

Yield  Stresses,  Observed  Yielding 

a 

none 

1 

^2 

14.5  ksl,  15.6  ksl 

Yield  Stresses,  .05%  criterion 

a  170  ksl 

1 

^2  15.6  ksl,  16,8  ksl 


REPORTED 
[Ref.  21 

31  msl 
20  msl 
6  msl 
.23 
.13 

216  ksl 
20  ksl 


49 


Table  3.  A  comprehensive  listing  of  B/Al  specimens  cut  from  panel  2 
as  shown  in  Figure  27. 


Specimen  and  Dimensions  in  inches  Notch 


Fiber 

Angle 

(t  X  W  X  1) 

Cut,  (in) 

Use 

16 

0® 

.0760 

X  .521  X  9.0 

— 

17 

30® 

.0780 

X  .502  X  10.0 

18 

45® 

.0780 

X  .527  X  9.0 

— 

19 

60® 

.0780 

X  .524  X  9.0 

— 

Elastoplastic 

20 

90® 

.0780 

X  .538  X  9.3 

— 

Characterization 

21 

0® 

.0775 

X  .500  X  9.0 

— 

22 

10® 

.0770 

X  .443  X  9.0 

— 

23 

20® 

.0780 

X  .465  X  8.0 

— 

24 

65® 

.0780 

X  .481  X  9.0 

— 

25 

45® 

.0780 

tapered  specimen 

— 

Preliminary 
plastic  zone 

test 

26 

90® 

.0775 

X  2.950  X  9.8 

.163 

Plastic  zone 
Measurements 

27 

60® 

.0775 

X  2.950  X  9.5 

.312 

at  60®  Plastic  zone 

Measurements 
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Table  A.  Mechanical  properties  of  B/Al  (tensile,  room  temperature) 
obtained  from  tests  outlined  in  Chapter  3  conducted  on 
specimens  17  through  24. 


Fiber 

Angle 

Elastic 

Modulus 

(msi) 

Elastic 

Poisson's 

Ratio 

Plastic 

Poisson's 

Ratio 

Ultimate 

Stress 

(ksi) 

Yield 

Stress 

(ksi) 

0“ 

27.05 

.232 

.577 

202.1 

none 

O 

o 

31.54 

.213 

.832 

70.4 

36.1 

20" 

26.96 

.278 

.757 

51.6 

19.7 

O 

O 

20.86 

.287 

.700 

34.2 

12.1 

45" 

15.99 

.260 

.646 

25.2 

12.7 

60" 

19.18 

.164 

.400 

17.6 

9.0 

65" 

20.77 

.656 

.393 

16.8 

11.7 

90" 

18.22 

.157 

.004 

22.2 

10.9 
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Table  5.  A  summary  of  the  failure  modes  and  ultimate  longitudinal 
failure  strains  for  specimens  17  through  24. 


Specimen 

Orientation 

Failure 

Mode 

Failure 

Strain 

0“ 

middle 

.0068 

o 

D 

end 

.0060 

20“ 

middle 

.0617 

30“ 

middle 

.0333 

45“ 

middle 

.0336 

60“ 

end 

.0078 

65“ 

end 

.0039 

90“ 

middle 

.0105 
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- ► 

FIBER  DIRECTION 


SIZE;  12"  X  24"  x  .077" 

MATRIX  /  FILAMENT:  6061  A1  /  5.6  B 
11  PLIES,  UNIDIRECTIONAL 
VOLUME  FRACTION:  47.5%  Fibers 


.  077’ 


Figure  1.  Specimen  cutting  plan  for  B/Al  panel  1.  A 

photograph  of  the  magnified  cross  section  shows  the 
11  ply  thickness.  The  section  was  cut  with  a 
diamond  saw  and  then  polished.  Note  that  several 
fibers  have  cracked  ends. 
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Figure  2.  Schematic  for  mechanical  properties  testing.  Load 
versus  strain  histories  were  recorded  on  X-Y 
plotters. 


LOAD 


5- 


STRESS,  (ksi) 


LONGITUDINAL  STRAIN,  (in/in) 


Figure  3.  Results  of  comparison  test  between  strain  gages  and 
the  clip  gage.  This  data  was  taken  from  B/Al 
specimen  5. 
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STRESS, 


Figure  4.  Typical  stress-strain  behavior  obtained  for  B/Al 
principal  material  directions  . 


STRESS,  (ksi) 


Figure  5.  Last  eight  load/unload  cycles  for  B/Al  specimen  7. 

Each  cycle  is  spaced  apart  for  clarity.  Note  the 
lack  of  residual  strain  when  unloaded.  The 
specimen  failure  is  indicated  by  the  x-mark 
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FSH)  ‘SS3HiS 


periodic  unloadings.  The  specimen  failure  is 
indicated  by  the  x-mark 


ELfiSTlC  MODULUS.  (MSI) 


22.00 


21  .GO 


20.0C-r  A 


19.00  + 


16.00  -f- 


^  A  ^ 


_L 


17.00  n - ^ - 1 - ^ - 1 

G.CC  2.0C  H.oo  6.00  E  .00 

PLRSTIC  STRAIN.  (X  .001 


10. CC  12. GO 


Figure  7.  Elastic  modulus  versus  plastic  strain  data  taken 
from  specimen  8.  The  modulus  is  observed  to 
decrease  substantially  with  Increasing  plastic 
strain.  See  Figure  6 
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Figure  8.  Schematic  for  fatigue  testing.  A  travelling 

microscope  was  used  to  observe  and  measure  fatigue 
crack  growth. 
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CRACK  INITIATION  IN  SPECIMEN  1 


0*  B/Al  3.0"  Wide  Notch  cut:  .168" 


DIRECTION 


♦ 


Figure  9.  Schematic  representation  of  cracks  in  B/Al 

emanating  from  the  notch  tip  and  propagating 
parallel  to  the  fibers. 


ber  angle  is  determined  from 
ge  after  cutting 
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NORMfiLIZEO  STRAINS 


BORON/ ALUMINUM.  CH  1  VALUES 


Figure  12.  Normalized  strain  versus  fiber  orientation  for  B/Al. 
The  normalized  shear  strain  peaks  around  40“ . 
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STRESS.  (PSD 


LONGITUDINAL  STRAIN 


Figure  14.  Stress-Strain  curves  for  a  hypothetical  material 
with  anisotropic  elastic  and  isotropic  plastic 
behavior,  generated  with  the  UNIAX.F  program. 


66 


TRANSVERSE  STRAIN 


Figure  15.  Transverse  versus  longitudinal  strains  for  a 

hypothetical  material  with  anisotropic  elastic  and 
isotropic  plastic  behavior,  generated  with  the  UNIAX.F 
program 
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STRESS.  (PSD 


LONGITUDINRL  STRRIN 


Figure  16.  Stress-strain  curves  for  a  hypothetical  material 
with  isotropic  elastic  and  anisotropic  plastic 
behavior,  generated  with  the  UNIAX.F  program. 


TRANSVERSE  STRAIN 


0.0000  .0100  .0200  .0300 

LONGITUDINRL  5TRRIN 


Figure  17.  Transverse  versus  longitudinal  strains  for  a 

hypothetical  material  with  isotropic  elastic  and 
anisotropic  plastic  behavior,  generated  with  the 
UNIAX.F  program. 
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FIBER  DIRECTION 


SIZE:  12"  X  24"  x  .077" 

MATRIX  /  FILAMENT:  6061  A I  /  3.8  B 
11  PLIES.  UNIDIRECTIONAL 
VOLUME  FRACTION:  47. 5X  Fibers 


Figure  18.  Specimen  cutting  plan  for  B/Al  panel  2. 
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COUPON  SPECIMENS 


THICKNESS  -  .077* 


on  of  a  typical  coupon  specimen  used  in 
outlined  in  Chapter  3 
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MTS  810 


Figure  20.  Schematic  of  the  test  setup  with  data  acquisition 
system. 
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STRESS.  (PSD 


1/2  INCH  WIDE  ALUMINUM  SPECIMEN 


Figure  2 


LONGITUDINfiL  STRAIN 


.  Stress  versus  strain  data  for  clip  gage 
calibration  for  half-inch  wide  specimens. 
Calibration  was  conducted  on  an  aluminum  specimen 
with  known  material  properties. 


LONGITUDINRL  STRfilN.  CLIP 


Clip  gage  strain  versus  strain  gage  strain  data  from  Figure  21. 

This  plot  shows  the  linear  relation  between  the 

two  strains  even  in  the  post-yield  region.  The 

plot  deviates  from  the  expected  45*  line, 

indicating  that  a  calibration  factor  is  required. 


STRESS.  PS  I  (XiO 


Figure  23. 


computer  modeline  of  iines  indie 

J^ing  ot  the  experimental  dats 


fiber 

the 
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O.OGO  .2GG  .HGO  .600  .800 

LONGITUDINAL  STRAIN  (XIO"^] 


Figure  24.  Stress  versus  longitudinal  strain  for  various  fiber 
orientations.  (Sane  data  as  in  Figure  23  shovm  on 
an  expanded  scale.) 


TRANSVERSE  STRAIN  IXIO  -2) 


Figure  25.  Transverse  versus  longitudinal  strain  for  various 
fiber  orientations.  The  solid  lines  indicate  the 
computer  modeling  of  the  experimental  data. 


.0100 


.0500 


c.occo 


_ L— _ i -  \ 

~1 - ^ ^ — 

.0200  .G3G0  .ChOC 


LCNGITUDINRL  STRRIN 


Figure  26.  Transverse  versus  longitudinal  strain  for  various 
fiber  orientations.  (Same  data  as  in  Figure  25 
shown  on  an  expanded  scale.) 


POISSONS  RATIO 


Figure  27.  Elastic  Poisson's  ratio  data  versus  fiber  angle 
solid  line  is  the  curve-fit  for  determining  the 
values . 


ELASTIC  MODULUS 


PLASTIC  POISSONS  RATIO 


Figure  29.  Plastic  Poisson's  ratio  versus  fiber  angle.  The 
dashed  line  Is  for  the  least-squares  curve-fit. 
The  solid  line  is  the  preferred  curve-fit. 
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EFFECTIVE  STRESS-STRAIN  FROM  LEAST-SQUARES  CURVE-FIT 


EFFECTIVE  PLHSfiC  STRfilN 


Figure  30.  Effective  stress  versus  effective  plastic  strain  for 
the  r-values  found  from  the  least-squares  curve- 
fit.  Note  that  the  0®  data  is  not  modeled  well. 
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EFFECTIVE  STRESS.  PSI 


Figure  31.  Effective  stress  versus  effective  plastic  strain  for 
the  preferred  curve-fit  r-values.  The  solid  line 
in  the  modeling  curve  . 


tapered  specimen 


Figure  32.  Schematic  of  the  tapered  specimen  with  dimensions 
and  strain  gage  locations.  The  specimen  was  cut 
using  electric  discharge  machining. 
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LOm.  LBS 


NET  SECTION  STRESS 


Figure  34.  The  data  presented  in  Figure  33  is  presented  here 

as  net  section  stress  versus  strain.  The  data  falls  nearly 
along  the  same  line.  Data  from  the  45*  coupon 
specimen  and  finite  element  modeling  is  included 
for  comparison. 
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Y 


REMOTE 

GAGE 


Figure  35.  Schematic  of  specimen  26  for  the  transverse  crack 
study.  The  strain  gage  strip  was  cut  into  two 
sections  and  positioned  as  shown. 
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Figure  37.  Schematic  of  specimen  27  for  the  oblique  crack 
study.  The  strain  gage  strip  was  cut  into  two 
sections  und  positioned  as  shown 
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rore%Ue%rrlng  the  1.5  overloaa. 


FIBERS 


0  *  GAGE  LOCATION  ANGLE 
r  «  GAGE  DISTANCE  FROM  CRACK 


Figure  39,  Global  coordinates  convention  for  the  study  of 
cracks  in  Chapter  4, 


Figure  40.  Plastic  zone  shape  at  the  crack  tip  generated  with 
ELLPLAS.F  program.  Fibers  are  perpendicular  to 
the  crack,  and  loading  is  as  shown.  The  crack 
length  is  .31"  (crack  shown  not  to  scale  with  plastic 
zones) 
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Figure  41.  Plastic  zone  shape  at  the  crack  tip  generated  with 
ELLPLAS.F  program.  The  fibers  are  parallel  to  the 
crack,  and  loading  is  as  shown.  The  crack  length 
is  .31"  (crack  shown  not  to  scale  with  plastic 
zones) 
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Figure  42.  Plastic  zone  shape  at  the  crack  tip  generated  with 
ELLPLAS.F  program  for  an  isotropic  material.  The 
crack  length  is  .31"  (crack  shown  not  to  scale 
with  plastic  zones). 
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SCALE:  2x  LIFE  SIZE 
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Figure  A3.  Plastic  zone  shape  at  the  crack  tip  generated  with 
ELLPLAS.F  program  for  fibers  parallel  to  the  crack 
and  load  angle  at  90*.  This  figure  is  similar  In 
nature  to  Figure  Al.  The  crack  length  is  .39" 
(crack  shown  no  to  scale  with  plastic  zones) 


O’  *  6560  psi> 
o  *  9840  psi- 
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Figure  AA.  Plastic  zone  shape  at  the  crack  tip  generated  with 
ELLPLAS.F  program.  The  fibers  are  parallel  to  the 
crack  with  the  load  angle  at  120*.  The  crack 
length  Is  .39"  (crack  shown  not  to  scale  with  plastic 
zones) .  This  crack  geometry  and  loading  is  similar 
to  that  of  specimen  27 


6560  psi 
9840  psi 


Figure  45.  Plastic  zone  shape  at  the  crack  tip  generated  with 
ELLPLAS.F  program.  The  fibers  are  parallel  to  the 
crack  with  the  load  angle  at  150®.  The  crack 
length  is  .39"  (crack  shown  not  to  scale  with 
plastic  zones) . 
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Figure  46.  Effective  stress  along  the  X-axis  for  specimen  26. 

Note  that  the  plastic  zone  boundary  is  determined 
by  where  the  plot  crosses  the  effective  yield 
line.  The  gage  locations  along  the  axis  are  as 
shown . 
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Figure  47.  Effective  stress  along  the  Y-axls  for  specimen  26. 

Jote  that  the  plastic  zone  boundary  Is  determined 
by  where  the  plot  crosses  the  effective  yield 
The  gage  locations  along  the  axis  are  as 

shown . 
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Figure  48.  Effective  stress  along  the  X-axis  for  specimen  27 
Note  that  the  plastic  zone  boundary  is  determined 
by  where  the  plot  crosses  the  effective  yield 
line.  The  gage  locations  along  the  axis  are  as 
shown . 


Figure  49.  Effective  stress  along  the  Y-axis  for  specimen  27. 

Note  that  the  plastic  zone  boundary  is  determined 
by  where  the  plot  crosses  the  effective  yield 
line.  The  gage  locations  along  the  axis  are  as 
shown.. 
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